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Abstract 

We propose a parametrization for the generahzed parton distributions (GPDs) 
which is based on representation of parton distributions as an infinite series of t- 
channel exchanges. The entire generahzed parton distribution is given as an infinite 
sum over contributions of generahzed hght-cone distribution amphtudes in the t- 
channel. We also discuss the relations of the lowest Mellin moments of GPDs to 
basic mechanical characteristics of the nucleon as a compound system. 



Introduction 

After the first experimental data on deeply virtual Compton scattering (DVCS) [jl], 0, |^, ^ 
and hard exclusive meson production |^ have been presented, it is the right time to ad- 
dress the problem of the extraction of the twist-2 generalized parton distributions (GPDs) 
from the observables. One of the possibilities to do this is to introduce a parametrization of 
GPDs and then to extract the corresponding parameters from the data. The parametriza- 
tion should respect all general constrains for GPDs (for reviews see |, ^, |1^) and it 
should have clear physical interpretation. 

To our best knowledge in all calculations of observables for hard exclusive processes 



the parametrization for GPDs based on double distributions (DDs) [TT|, |T2|, |T3|] has been 
used. In this parametrization (supplemented by the contribution of the D-term |14|) 
an important property of polynomiality of the Mellin moments of GPDs is satisfied 
automatically. The double distributions are modeled usually by the ansatz of the form 
|T3[| (for details see next section): 



F{(3,a) = h{a,(3)q{(3). (1) 

Here g(/5) is usual quark distribution and the profile function h{a,P) popularly is chosen 
usually in the form [|l3i: 



hia 3) = r(26 + 2) [jl-mr-a^t 
22^'+ir2(6+i) (i-\(3\)2b+i ■ 
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Here the parameter b characterizes the strength of the ^ dependence of the resulting GPDs. 
The advantage of the ansatz (|l]) is its simphcity, however this ansatz is too restrictive to 
have enough flexibihty in modehng of GPDs. Also this form of DD does not commute 
with the QCD evolution, i.e. if we assume this ansatz at one normalization point than at 
higher normalization points it is generically impossible to put the resulting DD into the 
form given by eq. |16 . 



In the present paper we suggest alternative way to parametrize GPDs based on the 
partial wave expansion of the GPDs. Before the discussion of the new parametrization 
we present in the next two sections the discussion of important property of polynomiality 
of Mellin moments and the relations of the lowest Mellin moments to the mechanical 
properties of the nucleon. 

In this notes we shall restrict ourselves only to the singlet quark GPDs. The formulae 
relevant to the nonsinglet GPDs are collected in Appendix A. We also do not consider 
the gluon GPDs because their consideration is very close in spirit to that of the quark 
singlet GPDs. For GPDs we use notations of X. Ji (see e.g. 0). Generically GPDs 
are functions of three variable (the fourth is log-dependence on the scale) H{x,^,t), to 
simplify notations we shall omit the variable t in the case when a quantity is assumed to 
be at t = or when the t-dependence is irrelevant for the discussion. 



Digression about polynomiality of Mellin moments and D-term 



In this section we do not write explicitly the t dependence of the GPDs as irrelevant 
for the present discussion but we always assume it. One of the non-trivial properties of 
the generalized parton distributions is the polynomiality of their Mellin moments. The 
polynomiality property means that 



jjxx^ E{x,i) = + ef ) + . . . + eS^;, 



(3) 



Due to the fact that the nucleon has spin 1/2, the coefficients in front of the highest power 



of ^ for the functions H and E are related to each other [15 



-h 



(TV) 
N+l 



(4) 



The polynomiality conditions (^ strongly restrict the class of functions of two variables 
H{x,C,) and E{x,C,)- For example the conditions (|^) imply that GPDs should satisfy the 
following integral constrains M: 



1 dx 

-1 X 



H{x,i + xz)-H{x,i) 



1 dx 
-1 X L 



E{x,i + xz)-E{x,i) 



E ^Si ^" • (5) 



n=0 



Note that the skewedness parameter ^ enters the Ihs of this equation, whereas the rhs 
of the equation is ^-independent. Therefore this ^-independence of the above integrals 
is a criterion of whether functions H{x,C,), E{x,C,) satisfy the polynomiality conditions 



2 



Simultaneously these integrals are generating functions for the highest coefficients 
npf^i- Also the condition shows that there are nontrivial functional relations between 
functions H{x,^) and E{x,^), more on this see below. 

An elegant possibility to implement the polynomiality conditions (^) for the GPDs is 
to use the double distributions [[l^, 0, 0. In this case the generalized distributions are 
obtained as a one-dimensional section of the two-variable double distributions F{j3, a) 
and K{P, a): 



H{x,0=f^d/3 da5{x- (3-aS) F{(3,a) , (6) 

J-l "'-1+1/31 



E{x,i)= d(3 da5{x- [3-ai) K{l3,a) . (7) 

J-l "'-1+1/31 



"1 rl 

-1+ 

and an analogous formula for the GPD E{x,C,)'- 

"1 rl 

-1+1/31 



[Recently in refs. [|T^, Q the inversion formula has been derived] 

It is easy to check that the GPDs obtained by reduction from the double distributions 
satisfy the polynomiality conditions (H) but always with /ij^i = ej^^^ = 0, i.e. the 
highest power of ^ is absent. In other words the parametrization of GPDs in terms of 
double distributions is not complete. It can be completed adding the so-called D-term to 
eq. ® 0: 



E{x,0 



df3 



dp 



1-1^ 

-1+1, 
1-1/^ 



1+1/31 



da 6{x- (3- ^iP, a) + 6 
da 6{x-(3-aO K{(3,a)-e 



X 



1 - 



X 



D 



X 



(9) 



Here D{z) is odd function having a support —l<z<l. In the Mellin moments the 
D-term generates the highest power of ^: 



uiN) 
"-N+l 



-e 



(N) 
N+1 



dz z^ D(z) . 



(10) 



Note that for both GPDs H{x,C,) and E{x,C,) the D-term is the same, it contributes to 
both functions with opposite signs. 

The D-term evolves with the change of the renormalization scale according to ERBL 



19, G0[ evolution equation. Hence it is useful to decompose the D-term in Gegenbauer 



series (eigenfunctions of the LO ERBL evolution equation): 

d, Cf/'(z) + 4 Cl^\z) + d, Cl^''(z) + 



D{z) 



l-z- 



d^ + d'iCl'\z)+d'^Cl'\z) + 



(11) 

(12) 
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GPDs and "mechanical properties" of the nucleon 

The x-moments of the GPDs H and E play a special role as they are related to the form 
factors of the symmetric energy momentum tensor. The nucleon matrix element of the 
traceless part of the symmetric energy momentum tensor is characterized by three scalar 
form factors [0, 15 1. These form factors describe the mechanical structure of the spin— | 



system-nucleon. At zero momentum transfer these three form factors can be related to 
three basic static mechanical characteristics of the composite system - nucleon. Two of 
them are energy-momentum and angular momentum carried by partons. The third one 
is related to the constants (ii(0) and (if (0) in expansion of the D-term, see eqs. ([TT| , p!2|) . 



These constants can be interpreted as the characteristic of the spatial distribution of the 
stress tensor or as of the spatial distribution of "forces" experienced by various species of 
partons inside of the nucleon. At t = we can express the x-moments of the GPDs as 



follows 15 : 



dx X {H{x, + E{x, 0) = 2 J« , (13) 

dx X H{x, = M^ + - di . 
J-i 5 

Here J*^ is a fraction of the nucleon angular momentum carried by all quarks, di is the 
first Gegenbauer coefficient in the expansion of the D-term (|lT]) and is a momentum 
fraction carried by quarks and antiquarks in the nucleon: 

= y dxx {q{x) + q{x)) . (14) 
q JO 

Analogous expressions one can write for gluon GPDs. To see physics interpretation of 
nucleon "mechanical" characteristics J^,M2 and di let us write their relations to the 
static energy-momentum tensor of the nucleon. In the Breit frame (in which the energy 
transfer A'' = 0) and non-relativistic limit (A^ -C mff) we can define static energy- 
momentum tensor of the nucleon as (see e.g. [^ ): 



1 r d'^A 

T^Ar) = ^ y ^""^ (P^ S\T^M\P, S) , (15) 

where T^(0) is the QCD operator of the symmetric energy momentum tensor of quarksQ 
and 5* is the polarization vector of the nucleon. Now it is easy to express the constants 
J^, M2 and di in the following way: 

= — !dhTg{f) 

= J e''^ Si r, Tg{r) . 
*An analogous equation and the discussion below can be trivially extended for gluons as well 
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These relations illustrate the interpretation of and J*^ as fractions of energy-momentum 
and angular momentum of the nucleon carried by quarks and antiquarks. This follows 
from the fact that T^{r) and Tq^^t) can be interpreted as a spatial distribution of energy 
and momentum carried by quarks correspondingly. 

The constant di can be expressed in terms of the stress tensor Tij{r): 



d, = J d\ T^{f) (^rV - i 5^^r2) . (16) 

If one would consider the nucleon as a continuous medium than T^j{r) would characterize 
the force experienced by quarks in an infinitesimal volume at distance r from the centre 
of the nucleon. More detailed distribution of forces one can obtain from the t-dependence 
of di{t) because of relation 




Obviously the sums + = 1, J*^-|-J*^ = l/2 and di + d'^ = d are scale indepen- 
dent. Two first sums corresponds to the total momentum and total angular momentum of 
the nucleon. The third constant d using eq. ([T6|) for the conserved total energy-momentum 
tensor can be represented in the form: 



rf=^ fd^rr'pir), (18) 
9 J 

where we introduced the following parametrization of the total (quarks+gluons) static 
stress tensor: 



sir 



(19) 



The function s{r) and p{r) are related to each other by conservation of the total energy- 
momentum tensor. The function p{r) can be interpreted as the radial distribution of the 
"pressure" inside the nucleon. We note also that the t-dependence of the GPDs provides 
us with more detailed spatial images of the nucleon, see e.g. |^ 0, ^ . 



The estimate which is based on the calculation of GPDs in the chiral quark soliton 
model at a low normalization point /x ~ 0.6 GeV, gives rather large and negative 
value of di ~ —4.0. The negative values of this constant has a deep relation to the 
spontaneous breaking of the chiral symmetry in QCD, see EO, IHO, P[. 



Partial wave decomposition of the GPDs 

Here we discuss the decomposition of GPDs in t-channel partial waves. Such kind of de- 
composition is useful for understanding of physical mechanism contributing to generalized 
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parton distributions. In a sense we attempt to model GPDs not by "deforming" ( "skew- 
ing") the forward parton distributions but rather representing the GPDs as an infinite 
sum of the distribution amphtudes in the t-channeL Using the analytical continuation of 
the t-channel exchange representation we shall try to relate two "dual" approaches to the 
GPDs. 

The partial wave decomposition in the t-channel for singlet GPDs H{x,C,,t) and 
E{x,^,t) can be written as the following formal series pB| : 



OO 71+1 



H{x,^,t) = j:j:Bruit)9il-^ (1-^ C^i^ PAt\ (20) 



n=l i=0 

odd even 



n=l !=0 

odd even 




OO 71+1 / 

E{x,u) = j:T.Cni{t)e d/U^] p,{^] . (21) 



Here the index / corresponds to an exchange in the t-channel with the orbital momentum /, 
because in the hard regime the t-channel scattering angle 9t is related to ^ via cos 9t = 1/^. 
The coefficients Bni{t) and Cni{t) are generalized form factors. Note that for the fixed 
power of the Gegenbauer polynomial n the orbital momentum can reach maximally the 
value / = n + 1, this is the consequence of the Lorentz invar iance [^. The formal series 
(pop corresponds to analytical continuation (crossing) of the corresponding expansion for 
the generalized distribution amplitudes |^ entering the description of hard processes like 
7*7 hh, see discussion in refs. [|26|, pTSl . 

Below we shall discuss mostly the GPD H{x,^,t) the expressions for E{x,^,t) are 
obtained by a trivial generalization. The form factors -B„„+i(t = 0) at zero momentum 



transfer are fixed in terms of the Mellin moments of usual singlet quark distributions ||26|| : 



Bnn+i{t = 0) = ——— / dx {q{x) + q{x)) . (22) 



2(n + 2) Jo 

The other coefficients Bni with I < n — 1 are not related to the forward distribution and 
thus are "truly non-forward" quantities. 

Each term in the sum (|20|) has a support — ^ < x < ^, however this does not imply that 
the GPD H{x,C,,t) has the same support because generically the sum ( PP] ) is divergent 



at fixed X and ^. In the formal solution ( pOD each individual term of the series gives the 
contribution to the amplitude which behaves as 1/^' at small ^ what would imply the 
violation of unitarity for the hard exclusive reactions. This "disaster" is avoided through 
the fact that the series is divergent and one can not calculate the asymptotic behaviour 
term by term. This situation is similar to duality idea: the entire generalized and usual 
quark distributions are given by infinite sum over t-channel exchanges. 

The formal representation (^) can be equivalently rewritten as the convergent seriesp]: 

OO 

H{x, e, t) = (1 - x') Y: ^n(e, t) C'J'ix) , (23) 

n=l 

odd 



t Actually this formula can be effectively used for numerical construction of GPDs from the coefficients 

Bnl 
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where 



2n + 3 



{n + l){n + 2) ^ 



(p+l)(p + 2) 



p+i 



odd 



!=0 
even 



Here Rnp{^) are polynomials in ^ of the order n introduced in ref. PSI, 53: 



p+n ^ 



r(^ + 



r(^ + i)r(|+p) 



n p 5 2 
2 + 2'2+^'^ 



(25) 



In the series (^) at small ^ the terms with I = n + 1 are dominant, so that in the limit 
^ — s> the forward distribution is recovered from eq. (pS]). Corrections to the forward 
limit of order at fixed x is completely determined by orbital momenta of / = n + 1 and 
I = n — 1. We see that the smaller skewedness ^ the less important the deviations of the 
the orbital momentum / from its maximal value of n + 1. We shall use this fact in next 
section to perform summation over orbital momenta. 

From either (|20|) or (|23|) we can easily obtain the partial wave decomposition of the 
Mellin moments 



(26) 



/•I ^"+1 /i\ T (^)t(N + l)(n + l)(n + 2) 
/ dx H{x, e, t) = V V BUt) Pi(-] [ , 

odd even V ^ / V ^ ^/ 



which obviously satisfies the polynomiality condition (^). For the lowest = 1 moment 
we have: 



5i2(t)--(5i2(t)-2i?io(t)) e 



j dx X H{x, ^, t) 
J dx X E{x, ^, t) 

which with help of relations ( ^2]) and (0) gives at t = 0: 

/I 1 / 1 2 \ 

_^ X H{x, = - - [^M^ - - fiio(O) j , 

f^dxxE{x,0 = 2J«-M2«-i(2J«-M2«-^Cio(0)) 



(27) 
(28) 

(29) 
(30) 



Prom this expression we can easily obtain the value of the first Gegenbauer coefficient of 
the D-term di at t = 0: 



di = -—M^ + B,o{0) 



(31) 



Equivalently we can extract the di using the fact that the D-term is the same (up to the 
sign) for H and E: 



rfi = -(2J«-M2'?)-Cio(0) 
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(32) 



Comparing eq. ( |3T[ ) with (^) we obtain new representation for the Ji's sum rule in terms 
of S-wave exchanges: 

J'^ = l (i?io(0) + Cio(O)) , (33) 
5 

which should be contrasted with the original formulation of this sum rule done in terms 
of D-wave exchanges: 

J'5 = |(5i2(0) + Ci2(0)) . (34) 
5 

We see that the can be extracted from the S-wave exchange in the t-channel. This is not 
surprising because between functions H{x,^) and E{x,C,) there are nontrivial functional 
relations, see e.g. eq. (|^). Generically for the A^-th Gegenbauer coefficient of the D-term 
we have two equivalent representations 

d^{t) = j:Bm{t) \ , (35) 



1=0 
even 



r(| + i)r(i) 



N+1 



:-irr(i + i 



1=0 
even 



This implies the following relations between coefficients BNi{t) and CNi(t): 

N+1 (_i)V2r(^i + i) N+1 (_i)V2r(i + i) 

V BNl{t) J- r T-T- = — ^ Cjyiit) -r- r r-v- . (37) 

r(i + i)r(i) h r(i + i)r(i) 



1=0 

even ^ ~ ^ ^" even 



Such kind of relations is useful for modeling of GPDs. 

The expression of the J*^ in terms of S-wave t-channel exchanges ( PHj) is useful because 
it allows us to analyze using various models and methods of low energy physics. 

Summing up the partial waves 

The partial wave decomposition of GPDs discussed in previous section can be used to 
write a parametrization of the GPDs in terms of "forward like" functions. To do this we 
introduce a set of functions whose Mellin moments generate the coefficients B^iit) and 

Cnl{t) 



Bnn+i-kit) = t dx Qk{x, t) , (38) 
JO 



Cnn+i-k{t) = dx x"^ Rk{x,t) , (39) 







where index k (always even number!) characterizes the deviation of the orbital momentum 
from the maximally possible value of / = n + 1. Note that the LO evolution of the 
functions Qk{x,t) and Rk{x,t) is governed by the DGLAP evolution equation. Since 
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-B„„+i(0) is fixed by the Mellin moments of forward distributions, see eq. (p2D, the function 
Qo{x,t = 0) is related to forward distributions: 



Qo{x) 



{q{x) + q{x)) - ^ f {q{z) + q{z)) 



(40) 



We sequence of functionsf] QQ{x^t),Q2{x,t),Q^{x,t), ... is introduced in such a way 
that the higher functions (with high index k) are more and more suppressed for small 
skewedness parameter ^. 



Using the methods developed in refs. [0, ^ we can derive the following integral 
transformation relating GPDs if(x, and functions Qk{y,t): 



CO ( 



H{x,^,t) = [H('\x,C,t)-H^'^\~x,C,t) 

fc=0 t ^ 
even 



(41) 



fc-3 



3/2 (x\ pJA f\....k-l-2 



odd 



where the function if*^'^^(x, ^, t) is given as the following integral transformation: 



1 -y^j Qkiy,t) 



H(>'\x,U) = 0{x>O - [ - 

TT Jyo y 



X 



1-k 



x1 — 2Xs — 

l-k 



x: 



.(42) 



Here Xs = '^j^z'^ ^-^id integration limits si,S2,S3 and y^ are given by the following 



expressions: 



= — : 



^2 = — : 



S3 



yo = 7^ 



1 
1 

1 
1 



i-ji-^2_2ii-xy)i-ji-e-eii-y 



i-ji-e + \2ii-xy)ii-ji-e~eii-y 



i + ji-e- W2(i-xy) 1 + vi-e -e(i-y 



X 1- vi-e 



2 i-ji-e 



(43) 



Actually the integral transformation (^) can be written compactly as: 
^The discussion below is applied to the functions Rk{x, t) as well 

^Remember that here we discuss only the singlet (C — +1) GPDs, collection of formulas for nonsinglet 
(C = —1) GPDs is given in the Appendix A 
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vr Jo y 



d \ 



ds 



X 



l-k 



xt 



2xs - e 



e{xl - 2x, - e) ■ (44) 



Derivation of eqs. (^,0) is given in the Appendix B. Also we note that the numerical 
realization of these equations is stable and fast, it takes fraction of a second on a PC to 
compute GPDs for given function Qk{x,t). 

The simple analytical application of these integral transformations is to estimate the 
small— ^, Regge-type, behaviour of the imaginary part of the leading order amplitude 
proportional to Y,kH^^\^,^). Taking the power-like parametrization for the function 
Qk{x) in the region of small y-Qk{y) ~ y'"", we get (for a > k — 1/2)^ 



k-c 



r(l/2)r(a-fc + l/2) 
r(a-A; + l) ' 



(45) 



which is in an agreement with the result of Ref. |3T|, ^ for A; = 0. Also we see that for 
the case when the a does not increase strongly with increasing of k the leading term is 
determined by the functions Qo{x) which is completely fixed by the forward distributions. 

The original idea of Refs. ^ to relate the GPDs to functions which are evolved 
according to usual DGLAP equation, as pointed out in Refs. |3^, implies nontrivial 
(^ dependent) constrains for the support of the "effective forward functions" . This makes 
the parametrization of these functions almost impractical. Our construction, given by 
eq. (|4^) , seems to be free of this problem. 

The leading order amplitude of hard exclusive reactions is expressed in terms of the 
following elementary amplitudeHJ]: 



dx H{x, C,, t) 



1 



+ 



1 



X 



^ + iO x + C.-iO 



(46) 



Now we can express the amplitudes in terms of "forward-like" functions Qk{x). For this 
we substitute the formal series representation for GPDs (|20| ) into eq. (^) and obtain the 
partial wave decomposition of the amplitude in the t-channel: 



OO 71+1 / 1 \ 

n=l 1=0 V^/ 



(47) 



odd even 



Substituting into this equation the expression for the coefficients Bni{t) in terms of func- 
tions Qk{x,t) (|39D, we can sum up the partial waves with the result: 



Jo X 



1-f + x^ 



+ 



2x 



2 S 



+ x' 



kO 



(4J 



^For the case of a < /c — 1/2 we obtain that i?'^'^^ ^) ~ with the coefficient which depends also 
on behaviour of the function Qk{y) at large y. 

II We restrict ourselves to the singlet or even signature amplitudes only. For the nonsinglet (odd 
signature) amplitude see Appendix A. 
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Using this relation we can write down explicit expressions for the real and imaginary parts 
of the amplitudes: 



dx °° 
e fc=o 



dx °° 

^ k=0 



dx °° 
— « — k=0 



^ - X2 - 1 



l~f + X^ 



+ 



1 + f + x^ 



2 5, 



kO 



1 + ^ + X2 



2(5, 



kO 



(49) 



These expressions would allow us to analyze the experimental data on hard exclusive 
reactions in terms of "forward-like" functions Qk{x, t) and Rk{x, t). Among these functions 
only Qo{x) is fixed in terms of forward distributions, see eq. (^). Note also that the 
smaller ^ the more suppressed the contribution of the functions Qk{x, t) and Rk{x, t) with 
higher k to the amplitude. We also note that from point of view of practical applications 
the numerical calculations of the amplitude according our formulae (|49D is more stable 
and faster than calculations of more singular integral (|46|). From expression (^91) for 
the imaginary part of the amplitude (corresponding to —7rH{^,^,t)) we see that for its 
calculation we need to know the functions Qk{x) (and hence forward parton distributions) 
for X > {1 — a/1 — C^)/^ > ^/2. In this way the problem (discussed in refs. [^) that 
in the parametrization of GPDs in terms of double distributions one needs to sample 
forward distributions down to very small (unmeasured) values of x is solved. In our 
parametrization we see explicitly that the contribution to the amplitude of very small x 
is strongly suppressed as it should be from physics point of view. 

From the expressions (^) one can easily obtain the behaviour of the amplitude at small 
^. Indeed, let us assume that at small x we have Qk{x,t) ~ a;""'-', then the contribution 
of each individual function Qk{x,t) to the amplitude has the form: 



ReA^'^) (e, t) ~ ImA('=) (e, t) tan |^![K_izLl)^ . (51) 

The first equation coincides with the small ^ behaviour of —tcH^''\^,^) (see eq. (|i5|)) as 
it should be. The second equation reproduces the result of general dispersion relations 
for the even signature scattering amplitudes at high energies Also we see that if ak 



does not grow drastically with k, the leading contribution comes from the function Qo{x) 
which is fixed in terms of the forward parton distributions. 

To derive eqs. (^,^) we note that at small values of the skewedness parameter ^ the 
expressions for the real and imaginary part of the amplitude (^) can be rewritten for the 
case ak — k > 1/2 as (after rescaling of the integration variables 2x/C, — > x): 
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Im A{U) 



fc=0 

oo 

E 

k=0 



dx 



k-Qk 



1 dx 

rj.a^-k+1 



1 



+ 



25, 



kO 



°° dx 



I j^Ok-k+l 



-25, 



kO 



X 



(52) 



(53) 



Performing integrals in the above expressions we obviously obtain eqs. ( |50| , |5TD . 

Using eqs. (|49D one can try to fit experimental data on hard exclusive processes adopt- 
ing a simple parametrization of the functions Qk{x,t) (for k > 2) borrowed from analysis 
of the DIS data: 



Qk{x, t) = Nk— (1 - x)^' (l + 7fcv^ + 6kx) , (54) 

and analogous form for the functions Rk{x). Where, in principle, all parameters ak,Pk, 
etc. are functions of the momentum transfer squared tf^ Introduced parameters for 
/c = 0, 2 at t = can be related to the basic mechanical characteristics of the nucleon- 
M^, and di by: 



dx X Qo{x) 



Ml 



2 1 



['dxxQ2{x) = d^ + ^M^, 
Jo 12 

pdxxRo{x) = -(2J'^-M^ 

f 

Jo 



dx X R2(x) 



(2 J<3 _ M^) - di 



(55) 



Note also that the leading order evolution of the functions Qkix,t) and Rk{x,t) is given 
by the usual DGLAP equations, so that for their evolution one can use fast and effective 
algorithms developed for the DIS, see e.g. Refs. ^ 

ak,Pk, etc. can be obtained from the positivity bounds for GPDs P7|, 



Further constrains on the parameters 

R31 derived in 



the most complete form recently by Pobylitsa |40 



"Minimal model 



Discussions in previous sections give us a hope that in the analysis of the data on hard 
exclusive processes we can restrict ourselves to the finite set of the "forward-like" functions 
Qk{x, t) and Rk{x, t). This set can be systematically extended with increasing the amount 

**The minimal possible choice of the t-dependence can be afc(t) — + a' t ^ where a' is the slope 
of the corresponding Regge trajectory. 
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of experimental data. The most interesting "mechanical" properties of the nucleon are 
contained in the lowest functions for k = 0,2. In this section we discuss possible minimal 
set of the "forward-like" functions. Clearly the choice with only A; = is too restrictive, 
in particular, such a choice would imply that H{x,^,t) = —E{x,^,t) which is too strong 
constrain. One may try to fit the data with functions Qo{x,t) and Q2{x,t) (and Ro{x,t), 
R2{x, t)) one of those is fixed by the forward parton distributions. In the model with only 
two types of generating functions Qo{x,t) and Q2{x,t) (and Ro{x,t), R2{x,t)) one can 
easily derive from the constrain ( pTj) the functional relation for the function R2{x,t): 

R2{x, t) = -Q2{x, t) + Qo{x, t) + R^{x, t)- — [Qo{z, t) + R^{z, t)] . (56) 

Jx Z 

We see that in such model one has only two new "forward-like" functions {Q2{Xit) and 
RQ{x,t)), QQ{x,t) is essentially fixed by the forward distributions, see eq. (|40|)P], and 
R2{x,t) by the relation (FBI). 



)5 



In this paper we do not present our numerical results for GPDs in the "minimal model 
just discussed, the corresponding results we shall present elsewhere. Our main aim here 
was to set the theoretical framework for modeling of GPDs. 



Conclusions and outlook 



We analyzed the representation of the generalized parton distributions (GPDs) as the 
infinite sum over distribution amplitudes in the t-channel. Such representation is close in 
its spirit to the duality idea-the entire parton distribution is obtained as an infinite sum 
over distribution amplitudes in the t-channel. On basis of this dual representation we sug- 
gested a parametrization of the GPDs which satisfy automatically all general constrains- 
forward limit, polynomiality, etc. In this parametrization the GPDs and the leading 
order amplitudes are expressed in terms of "forward-like" functions Qk{x,t) {Rk{x,t)). 
The dependence of these functions on the scale in the leading order is governed by the 
DGLAP evolution equation. The index k which enumerates the set of functions Qk{x,t) 
{Rk{x,t)) describes the deviation of the orbital momentum / in the t-channel from its 
maximal value of / = n + 1 at fixed conformal spin n of the twist-2 operator. The basic 
mechanical characteristics of the nucleon-momentum and angular momentum fractions 
carried by quarks and additionally the radial distribution of forces experienced by quarks 
in the nucleon-are contained in the lowest functions Qo{x,t),Q2{x,t) {Ro{x,t), R2{x,t)). 

With our parametrization of GPDs still there are several points which should be worked 
out further: 

• To prove that suggested parametrization covers all possible forms of GPDs and 
derive the inversion formula expressing functions Qk and Rk through GPDs H and 
E. 

• To extend the analysis to the next to leading order. On NLO evolution of GPDs 
see, e.g. ref. 

An analogous expression can be also written for Ra{x) where in rhs the role of forward distribution 
is played by E{x, ^ — 0). 
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Appendix A 

In the main text we discussed only the singlet (C = +1) GPDs which are reduced to 
[q{x) + q{x)] in the forward limit. In this Appendix we collect formulae for the nonsiglet 
{C = —1) GPDs which in the forward limit correspond to [q{x) — q{x)]. Below we simply 
list the corresponding formulae, skipping the discussion. 
Partial wave decomposition: 



even odd 

.(..^,, ^ f |c„,w<- (1-3(1 -P)c-(|jfl 

even odd 



2\ /.\ /^N 



(57) 
(58) 



Mellin moments : 

^"+1 /i\ T (^)t(N + l)(n + l)(n + 2) 

dx X- Hix,^,t) = E E BUt) P. 1 ) . /x^/J: J ■ (59) 



n = l 1 = 

even odd 



A J 2"r(^ + i)r( 



N+n I 5 
2 2 



The lowest x^ Mellin moments of GPDs are related to the electromagnetic form factors 
of the nucleon 



and 



f^dx H{x,U) = F,{t) = '^Boi{t) 
dx Fix, e, t) = F2{t) = ^Coi(t) : 



Generating functions Qk and Rk are introduced in the same way as in eqs. ( P^D only in 
the nonsinglet case the number n is even. The function QQ{x,t = 0) is related to the 
nonsinglet forward quark distribution: 



(g(x) - q{x)) 



X dz 



2 (Qiz) -q{z)) 



(60) 
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Integral transformations Qk — ^ H. 



k=0 
even 



(61) 



k-2 

1=0 
even 



3/. l-\ pjl 



a Jo 



dy y'-'-' Qk{y,t) 



with H^''^ given by eq. 
Amplitudes 

The amplitude with C = — 1 exchange in the t-channel (odd signature) defined as 



A{^,t) = / dxH{x,^,t) 
Jo 



1 



1 



X — ^ + iO X + ^ — iO 
can be expressed in terms of the nonsinglet functions Qk{x,t) as follows: 



(62) 



Re A(e,t) 



1 



2a; 



+ 



1 



X 



k=0 



1 + ^ + x' 



■' € ^ fc=0 



2x_^2_l 



(63) 



From these expressions we can derive analytically the results for the small ^ behaviour of 
the amplitude for the case of Qk{x) ~ x""*. The result is 



'^Nfc Y(i/2)T{ak-k + l/2) 



ReA^''\^,t) ~ ImA(^)(^,t) cotan 



T{ak-k+l) 
( 7r(afc -k-l) 



(64) 



The first equation again coincides with the small ^ behaviour of —ttH^'^^^,^) as it should 
be. The second equation reproduces the well-known result of general dispersion relations 
for the odd signature scattering amplitudes at high energies. 

Appendix B 

Here we sketch the derivation of the eqs. (|4l|^j4^). The basic relation we use is the 
following^: 

■'■■'■We define the discontinuity as disCz=x.f{z) ~ ^ [f{x — iO) — f{x + iO)]. 
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1 

y 



' dy) 1-1 



X 

e 



Here 



z.. = 2- 



z- 



with < y < 1. To prove eq. ( p5| ) we note that: 



(65) 



(66) 



■TV 



1 



r(iv) 



Xq 



with 



x.. = 2 



il-s^)y 



This simple relation gives us 



(67) 



(6^ 



disc. 



ds z\ 



-N 



9 1- 



X 



,N 



d 



N-l 



^2 j 2N^NT{N) \ds^ 



(69) 



--x/i 



The step function in the above equation indicates that the zero of the 5-function should 
be inside the interval — 1 < s < 1. Further with help of the identity: 



+ .;- = (iv + i).;-, 

and the Rodrigues formula for the Gegenbauer polynomials: 



(70) 



(1-.=,C-,.,^ (-1,-^(1 



N-l 



(1-X^) 



2\N 



we arrive to the eq. (pSl). 

Now we consider the following function: 



21-1/2 



(71) 



(72) 



and let us compute its discontinuity with help of eq. (^) and generating function for the 
Legendre polynomials: 



disc_F('=)(.,,) = (i - ^)^(i - ^) (f) (^) y'^'^' 



(73) 
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In this expression we immediately recognize the integral kernel for the formal expansion 



with 



Now the trick is that we can compute the discontinuity of the function F^''\z, y) given by 
eq. (F^) in a different way. Namely, we take contributions to the discontinuity from the cut 



1 — \/l — C,'^ < Zs < 1 + -y/l — C,'^ and from the poles at = for > 2. The calculations 
are straightforward, the cut contribution gives us the expression (H), whereas the pole 



contribution provides us with the second line in eq. pl|). Further, analyzing the solutions 
of the algebraic equation — 2xs + = we arrive at eq. (^2|) . 
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